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Abstract: The Pierre Auger Observatory allows the measurement of both longitudinal profiles and 
lateral particle distributions of high-energy showers. The former trace the overall shower development, 
mainly of the electromagnetic component close to the core where the latter reflect the particle densities 
in the tail of the shower far away from the core and are sensitive to both the muonic and electromagnetic 
components. Combining the two complementary measurements, predictions of air shower simulations are 
tested. In particular the muon component of the tank signals, which is sensitive to hadronic interactions 
at high energy, is studied with several independent methods. Implications for the simulation of hadronic 
interactions at ultra-high energy are discussed. 
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Introduction 

During the last decade, air shower simulation 
codes have reached high enough quality that there 
is good overall agreement between the predicted 
and experimentally observed shower characteris- 
tics. The largest remaining source of uncer- 
tainty of shower predictions stems from our lim- 
ited knowledge of hadronic interactions at high en- 
ergy. Hadronic multiparticle production has to be 
simulated at energies exceeding by far those ac- 
cessible at terrestrial accelerators and in regions of 
phase space not covered in collider experiments. 
Therefore it is not surprising that predictions for 
the number of muons or other observables, which 
are directly related to hadron production in show- 
ers, depend strongly on the adopted hadronic inter- 
action models [ 1 ] . 

In this work we will employ universality features 
of the longitudinal profile of the electromagnetic 
shower component to combine fluorescence detec- 
tor and surface array measurements of the Pierre 
Auger Observatory. Using the measured depth 
of shower maximum, X max , the muon density at 
ground is inferred without assumptions regarding 
the primary cosmic ray composition. This allows a 
direct test of the predictions of hadronic interaction 
models. 



Parameterisation of surface detector 
signal using universality 

Universality features of the longitudinal profile of 
showers have been studied by several authors [2]. 
Here we exploit shower universality features to 
predict the surface detector signal expected for 
Auger Cherenkov tanks due to the electromagnetic 
and muonic shower components at 1000 m from 
the shower core. In the following only a brief 
introduction to the method of parameterising the 
muonic and electromagnetic tank signals is given. 
A detailed description is given in [3]. 

A library of proton and iron showers covering the 
energy range from 10 to 10 20 eV and zenith an- 
gles between 0° and 70° was generated with COR- 
SIKA 6.5 [4] and the hadronic interaction models 
QGSJET 11.03 [5] and FLUKA [6]. For compari- 
son, a smaller set of showers was simulated with 
the combinations QGSJET II.03/GHEISHA [7] 
and SIBYLL 2.1/FLUKA [8, 9]. Seasonal models 
of the Malargue molecular atmosphere were used 
[10]. The detector response is calculated using 
look-up tables derived from a detailed GEANT4 
simulation [11]. 

Within the library of showers, the predicted sur- 
face detector signal for the electromagnetic com- 
ponent of a shower at the lateral distance of 1000 m 
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is found to depend mainly on the energy and the 
distance between the shower maximum and the 
ground (distance to ground, DG — A" groun d — 
I max ). Here the signal of electromagnetic shower 
component is defined as that of all shower parti- 
cles except muons and decay products of muons. 
The signal at 1000 m depends only slightly on the 
mass of the primary particle (13% difference be- 
tween proton and iron primaries) and the applied 
interaction model (~ 5%). The functional form, 
however, is universal. The situation is similar for 
the expected tank signal due to muons and their 
decay products. In this case the shower-to-shower 
fluctuations are larger and the difference between 
proton and iron showers amounts to 40%. 




gQGSii.p is the muon signal p re( jicted by QGSJET 
II for proton primaries. The relative importance 
of the electromagnetic and muonic detector signal 
contributions at different angles is shown in Fig. 1. 

Constant-intensity-cut method 

Within the current statistics, the arrival direction 
distribution of high-energy cosmic rays is found to 
be isotropic, allowing us to apply the constant in- 
tensity cut method to determine the muon signal 
contribution. Dividing the surface detector data 
into equal exposure bins, the number of showers 
with ,5(1000) greater than than a given threshold 
should be the same for each bin 
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Using the independently measured mean depth of 
shower maximum (X max ) [12] the only remain- 
ing free parameter in Eq. (2) is the relative num- 
ber of muons iV™ 1 . For a given energy E, Njf- 
is adjusted to obtain a flat distribution of events in 
sin 2 8. 



Figure 1 : Electromagnetic and muon contributions 
to the detector signal as a function of zenith an- 
gle. Results of QGSJET II/FLUKA simulations 
are shown for 10 19 eV showers. 

After accounting for geometrical effects such as 
the projected tank surface area, the proton-iron av- 
eraged electromagnetic shower signal is parame- 
terised in dependence on the energy E, distance to 
shower maximum DG, and zenith angle 8. The 
difference between proton and iron shower pro- 
files is included in the calculation of the systematic 
uncertainties later. Similarly the universal shape 
of muon signal profile is parameterized simulta- 
neously for all model primaries, taking the overall 
normalisation from proton showers simulated with 
QGSJET II/FLUKA. The expected detector signal 
at 1000 m can then be written as 

Smc(E, 8, Xmax) = S em (E,8, DG) 

+Nf S ? GSII ' P (10 19 eV, 8, DG),(l) 

where AT, is the number of muons relative to 
that of QGSJET proton showers at 10 19 eV and 
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Figure 2: Sensitivity of the constant-intensity-cut 
method to the muon number for E = 10 19 eV. 

The sensitivity of this method to the muon number 
parameter in Eq. (1) is illustrated in Fig. 2. The 
best description of the data above 10 19 eV requires 
N™ 1 = 1.63. However, this result was obtained 
by using the measured mean depth of shower max- 
imum [12] in Eq. (1). Shower-to-shower fluc- 
tuations in AT max and the reconstruction resolu- 
tion cannot be neglected and have been estimated 
with a Monte Carlo simulation. Accounting for 
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fluctuations and reconstruction effects, the rela- 
tive number of muons at 10 19 eV is found to be 
L45±0.11(stat)tg;j$(sys). 
Knowing the muon number and the measured 
mean depth of shower maximum, the signal size 
at 9 = 38° can be calculated 

5 38 (10 19 eV) = 37.5 ± 1.7(stat):±2i(sys) VEM. 

(3) 

This value of S38 is a measure of the energy scale 
of the surface detector which is independent of 
the fluorescence detector. It is within the system- 
atic uncertainties of the energy determination from 
fluorescence detector measurements, including the 
uncertainty of the fluorescence yield [13]. It corre- 
sponds to assigning showers a ~ 30% higher en- 
ergy than done in the fluorescence detector-based 
Auger shower reconstruction (E = 1.3-Efd). 

Hybrid event and inclined shower anal- 
ysis 

Hybrid events that trigger the surface detector ar- 
ray and the fluorescence telescopes separately are 
ideally suited to study the correlation between the 
depth of shower maximum and the muon density at 
1000 m. However, the number of events collected 
so far is much smaller. For each individual event 
the reconstructed fluorescence energy and depth of 
maximum are available and the expected 5(1000) 
due to the electromagnetic component can be cal- 
culated directly. The difference in the observed 
signal is attributed to the muon shower component 
and compared to the predicted muon signal. 

For this study, high-quality hybrid events were se- 
lected for which the shower maximum was in the 
field of view of a telescope, 9 < 60°, and the 
Mie scattering length was measured. Furthermore 
the distance between the telescope and the shower 
axis was required to be larger than 10 km and the 
Cherenkov light fraction was limited to less than 
50%. The surface detector event had to satisfy the 
T5 selection cuts which are also applied in [13]. 

In Fig. 3, we show the muon signal derived from 
these hybrid events as function of distance to 
ground. The relative number of muons at 10 19 eV 
is found to be 

<U =1 . 3EPD = 1-53 ±0.05 
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Figure 3: Reconstructed and predicted muon tank 
signal contribution in dependence on the distance 
to ground for vertical and inclined hybrid events. 
The muon profiles expected from QGSJET II sim- 
ulations are indicated by the red (proton showers) 
and blue (iron showers) points for the energy scale 
E = 1.3-Efd- 

Nf\ E=EFD = 1.97 ±0.06, (4) 

consistent with the analysis above. 

A similar study has been performed for inclined 
hybrid events (60° < 9 < 70°). Within the limited 
statistics, good agreement between muon numbers 
of the inclined and the vertical data sets is found, 
see Fig. 3. 

In Fig. 4 we compare the results of the differ- 
ent methods applied for inferring the muon den- 
sity at 1000 m from the shower core. The rel- 
ative number of muons is shown as function of 
the adopted energy scale with respect to the Auger 
fluorescence detector energy reconstruction. Only 
the constant-intensity-cut method is independent 
of the energy scale of the fluorescence detector. 
Very good agreement between the presented meth- 
ods is found. 

Discussion 

Assuming universality of the electromagnetic 
shower component at depths larger than X niax , we 
have determined the muon density and the energy 
scale with which the data of the Auger Observatory 
can be described self-consistently. The number of 
muons measured in data is about 1 .5 times bigger 
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Figure 4: Comparison of the results on the relative 
muon multiplicity at 10 19 eV from different meth- 
ods. 

than that predicted by QGSJET II for proton show- 
ers. Consistent results were obtained with several 
analysis methods. 

The QGSJET II and SIBYLL 2.1 predictions for 
iron showers correspond to relative muon numbers 
of 1.39 and 1.27, respectively. Therefore, inter- 
preted in terms of QGSJET II or SIBYLL 2.1, the 
derived muon density would correspond to a pri- 
mary cosmic ray composition heavier than iron, 
which is clearly at variance with the measured 
X max values. The discrepancy between air shower 
data and simulations reported here is qualitatively 
similar to the inconsistencies found in composition 
analyses of previous detectors, see, for example, 
[14, 15, 16]. 

Finally it should be mentioned that the results of 
this study depend not only on the predictions of 
the hadronic interaction models but also on the re- 
liability of the model used for calculating the elec- 
tromagnetic interactions (EGS4 in this study [17]). 
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